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I .  BACKGROUND  AND  OBJECTIVES 


Xn  July  1976,  a  contract  was  awarded  by  tha  Office  of  Naval 
Rasaarch  to  Elactric  Boat  Division  in  tha  araa  of  computer 
analysis  of  waldad  ship  structuras.  Tha  overall  obj active  of 
this  rasaarch  was  to  apply  and  furthar  davelop  whera  necessary 
tha  tachnology  for  pradicting  distortion  of  submarine  structure 
dua  to  welding.  The  underlying  motivation  in- providing  this 
technology  is  that  it  can  be  used  In  design  and  construction 
of  submarines  to  help  determine  tha  most  cost  effective  weld 
joint  designs  and  welding  procedures.  The  specific  objectives 
of  tha  initial  contract  ware; 

1)  implementation  of  MXT  finite  element  computer 
programs  for  prediction  of  temperatures,  residual 
stresses,  and  distortion  dua  to  welding, 

2)  evaluation  of  MIT  computer  programs  on  a  fusion 
weld  of  two  plates  and  comparison  with  results 
obtained  by  MARC-CDC  finite  element  program  (1)*, 

3)  analysis  of  a  butt  waldad  unrestrained  plate 
with  MXT  computer  codes  and  comparison  with 
test  results. 


^Numbers  in  parenthesis  indicate  references. 


The  results  of  that  study  are  documented  in  Reference  2.  In 
June  1977,  a  contract  modification  extended  the  scope  of  the 
current  contract  to  the  following  major  technical  objectives: 

1)  transient  thermoplastic  analysis  of  a  butt  welded 
unrestrained  plate  using  the  Electric  Boat 
Division  (EBDiv.)  modified  NONSAP  finite  element 
computer  code, 

2)  further  development  and  application  of  EBDiv. 
Shrinkage  Force  Method  (SFM)  to  welding  analysis 
including  investigation  of  weld  bead  sequencing, 
non-symmetry  in  the  weld  groove  and  remelt  effects. 
The  SFM  is  a  simplified  analytical  approach  for  pre¬ 
dicting  weld  distortion.  The  method  is  described  in 
Section-  IV- A  of  the  report. 


The  results  of  this  study  are  documented  in  this  report. 

The  technical  objectives  for  the  modified  contract  were  developed 
on  the  basis  of  the  results  of  the  first  year's  work  and  dis¬ 
cussions  with  cognizant  ONR  and  MIT  technical  personnel. 
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Several  problem  areas  In  the  use  of  the  MIT  computer  programs 
for  welding  simulation  were  uncovered  during  the  first  year. 

As  a  result  of  these  problems,  EBDiv.  began  modification  of 
its  own  in-house  software  for  analysing  weld  problems.  This 
resulted  in  the  implementation  of  two  computer  codes.  The 
first  code,  TEMP  (2),  performs  the  nonlinear  heat  transfer  of 
two  dimensional  solids.  TEMP  was  used  in  the  first  year  of 
the  contract  to  calculate  the  transient  temperature  solutions 

I 

for  the  butt  welded  unrestrained  plate  problem.  The  second 
code,  NONSAP  (2),  was  modified  to  handle  two  dimensional  thermo- 
plasticity  problems  with  temperature  dependent  material  properties. 

I 

The  first  objective  of  the  modified  contract  was  to  study  the 
distortion  and  stresses  resulting  from  the  butt  welding  of  the 
unrestrained  plate  using  the  modified  NONSAP  program.  The 
transient  temperature  solutions  from  the  TEMP  program  would  be 
used  as  the  driving  thermal  gradients  for  the  problem. 

It  became  apparent  in  the  early  stages  of  this  study  that  the 
high  level  of  technology  available  in  the  more  sophisticated 
i  computer  programs  for  welding  simulation  would  be  costly  to 

;  employ  on  many  shipyard  problems.  This  is  particularly  true 

if  parametric  studies  involving  several  variables  were  required. 

j  This  degree  of  sophistication  would  be -needed  for  detailed 
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analyses,  such  as  residual  stresses  and  fracture  of  weldments.  '  n 

However,  many  shipyard  problems  stem  from  the  geometric  distortions 
caused  by  welding;  in  many  cases  design  and/or  fabrication 
decisions  can  be  made  on  the  basis  of  relative  distortion 
comparisons  through  parametric  study. 

In  light  of  these  current  requirements,  a  second  objective, 
the  study  and  application  of  approximate  "weld  simulation" 
methods,  was  established.  In  particular  the  Shrinkage  Force 
Method  (SFM)  has  been  extended  to  two  dimensional  solids.  The 
method  is  employed  in  the  GBNSAM  (3)  computer  program  and 
studies  concerning  weld  sequencing,  weld  symmetry,  and  remelt 
have  been  performed  under  the  modified  contract.  Electric  Boat 
has  concurrent  with  this  contract,  initiated  an  I RAD  (4) 

(Internal  Research  and  Development)  task  which  emphasizes  the 
employment  of  approximate  'Veld  simulation"  techniques  such 
as  SFM  in  shell  welding.  Under  this  task,  the  weld  problem  is 
being  viewed  as  a  layered  shell  analysis  in  the  weld  region,  " 
with  the  layers  having  the  material  and  thermal  characteristics 
of  the  weld  deposit.  The  analytical  approach  utilizes  a  finite 
difference  shell  model  as  opposed  to  the  two  dimensional  solid 
finite  element  method  which  is  currently  employed  in  the  (  R 
studies. 


* 
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II .  EXPERIMENTAL  RESULTS  •  UNRESTRAINED  PLATES* 

II-A.  Test  Set-up 

Figures  II>1  and  II- l  contain  plan  and  elevation  views  of  the 
test  sat -up  along  with  a  detail  of  the  weld  groove  geometry. 

A  total  of  IS  weld  beads,  3  initial  root  passes  by  shielded 
metal  arc  and  12  fill  passes  by  submerged  arc  were  made.  Fig¬ 
ure  IV- 3  contains  a  macro  section  of  the  final  vrelded  region. 

Bead  number  one  was  removed  by  baekgrinding. 

Distortion  data  was  collected  via  dial  indicators,  potentiometer 
transducers  and  strain  gage  beam  instruments .  Their  relative 
locations  are  given  in  Figure  II- 3. 

I I -B ,  Distortion  Data  Analysis 

Figures  II-4A  and  4B  contain  the  measured  deflections  at  each 
of  the  gauge  locations  for  welding  of  the  first  side.  There 
is  a  significant  discrepancy  between  the  readings  from  the 
transducer  labelled  B  and  the  five  other  gauges.  Being  unable 
to  account  for  the  discrepancy,  it  was  decided  that  the  gauge  B 
data  would  not  be  used  in  the  evaluation  of  experimental  results. 

*This  test  was  performed  under  a  197S  EBDiv.  IRAD  project  (Ref.  S) . 


Figurt  II* S  shows  tha  average  daflactions  of  gaugas  A  and  D, 

C  and  F  as  compartd  with  the  transducar  at  E,  tha  early  baad 
distortion  response  at  E  is  graatar  than  tha  average  of  tha 
rasponsas  at  tha  ands  of  plata.  This  indicatas  a  slight  bowing 
upward  which  ravarsas  to  downward  at  about  baad  7.  Overall, 
Pigura  II* S  shows  a  naarly  uniform  daflaction  pattarn  along 
tha  langth. 

Pigura  11*6  contains  tha  plots  of  changes  of  daflactions  for 
aach  baad.  Tha  data  indicate  that  tha  largest  changes  occur 
in  going  from  baad  1  to  baad  3,  with  tha  second  largest  incre¬ 
ment  occurring  in  tha  deposit  of  baad  5. 

Figure  11*7  contains  plots  of  tip  daflactions  and  rotation  at 
E  along  with  their  associated  increments.  Tha  data  in  Fig- 
ura  11-7  was  usad  as  the  basis  fcT  comparison  with  analytical 
rasults. 
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III.  MATERIAL  PROPERTY  DESCRIPTION 


III-A.  Introduction 

One  of  the  most  difficult  aspects  of  attempting  to  analytically 
characterize  the  welding  process  is  the  definition  of  material 
properties  over  the  wide  range  of  temperatures  involved.  For 
many  existing  steels,  this  definition  does  not  appear  to  exist. 

An  additional  complication  occurs  when  the  material  properties 
of  interest  are  sensitive  to  both  temperature  and  the  rate  of 
cooling  (heating) .  Existing  available  data  is  very  limited  in 
scope,  even  for  the  most  recently  developed  high  strength  steels. 

A  successful  analytical  simulation  of  welding  is  directly 
dependent  on  an  ability  to  characterize,  through  a  mathematical 
material  model,  the  behavior  of  the  material  in  the  weld  region. 
This  mathematical  model  is  in  turn  dependent  upon  properties 
such  as  the  elastic  modulus  and  coefficient  of  expansion,  etc. 
which  must  be  experimentally  determined.  There  is  a  degree  of 
doubt  as  to  the  existence  of  an  adequate  mathematical  material 
model  that  would  precisely  represent  the  behavior  in  the  weld 
puddle  region. 

In  this  study,  the  material  is  characterized  in  two  ways.  For 
the  shrinkage  force  analysis,  the  material  is  assumed  to  be 


IS 


quasi-linear.  The  properties  take  on  different  discrete  values 
over  certain  temperature  ranges  and  behavior  is  assumed  to  be 
linear  over  each  temperature  range.  The  NONSAP  analysis  employs 
a  thermal  plasticity  model  which  allows: 

a)  temperature  dependence  of  Young's  Modulus  (E) , 

Yield  Stress  (by) ,  Strain  Hardening  Modulus  (ET) 
and  Coefficient  of  Thermal  Expansion  (o^.) , 

b)  nonisothermal  Von  Mises  yield  function  and  its 
associated  flow  rule, 

c)  Isotropic  hardening  or  elastic-perfectly  plastic 
behavior. 

Both  of  these  models  assume  no  dependence  of  material  properties 
on  cooling  rates. 


III-B.  Properties  of  HTS 


A  number  of  different  sources  were  examined  in  an  attempt  to 
obtain  a  description  of  the  material  properties.  These  sources 
included  ASTM  publications,  military  specifications,  major 
steel  companies,  universities,  government  laboratories  and 
private  corporations.  In  most  cases  the  chemical  composition 
and  heat  treatment  methods  had  to  be  compromised  to  find  a 
similar  material  with  available  data.  In  all  cases  the  data 
was  incomplete  for  the  temperature  range  of  interest.  However, 
the  survey  was  useful  in  that  several  organisations  that  could 
perform  the  proper  testing  were  identified. 

The  HTS  plate  used  in  the  experiment  falls  under  ASTM  desig¬ 
nation  S37  class  1.  Class  1  specifies  a  normalized  heat  treat¬ 
ment  as  opposed  to  quenched  and  tempered.  Normalizing  requires 
a  ferrous  alloy  to  be  heated  to  a  suitable  temperature  above 
the  transformation  range  and  then  cooled  in  air  to  a  temper¬ 
ature  substantially  below  the  transformation  range.  The 
transformation  range  is  defined  by  those  temperatures  within 
which  an  austenite  structure  forms  during  heating  and  transforms 
during  cooling.  The  rate  of  cooling  controls  the  final  grain 
structure  obtained.  The  transformation  range  for  HTS  was 
taken  as  510 ®C  to  "60*C  and  the  major  property  changes  that 
occur  are  reflected  in  the  variations  that  occur  in  this  range. 
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The  following  discussion  will  bo  broken  down  by  individual 
thermomechanical  proparty.  Those  proportios  at  best  roprosont 
a  compromise  from  a  number  of  sources. 


Elastic  Modulus 


Data  supplied  by  the  United  States  Steel  Corporation 
Research  Laboratory  for  steels  with  austenitic  grain 
structure  show  a  variation  of  E  at  room  temperature 
of  from  1.9306  x  IQ12  dynes/cm2  (28  x  106  psi)  to 
2.068S  x  1012  dynes/cm2  (30  x  106  psi).  Some  general 
data  on  the  E  variation  with  temperature  was  supplied 
by  USS  but  it  did  not  cover  the  full  temperature 
range.  The  data  presented  in  References  6,  7  for 
HT80  steel  was  felt  to  be  more  accurate  and  therefore 
was  employed  in  this  study.  Figure  III-l  contains 
the  assumed  variation  of  elastic  modulus  with  respect 
to  temperature  for  HTS. 


Data  supplied  in  Reference  S  for  for  carbon  steel 
was  used.  This  data  was  available  for  temperatures 
less  than  6SQ*C  (1200*F)  and  extrapolated  far  beyond 
this  temperature.  The  thermal  coefficient  was  assumed 
to  increase  linearly  from  the  last  known  data  points 
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until  760*C.  At  this  point  it  was  assumed  to  drop 
to  a  value  that  represented  501  of  the  total  change 
experienced  in  going  from  room  temperature  to  760 °C. 
Although  this  drop  is  characteristic  of  steels,  the 
amount  of  decrease  was  completely  arbitrary  since  no 
data  was  available.  This  was  simply  an  attempt  to 
realistically  reflect  a  known  fact  about  the  aT 
behavior.  From  760*C  to  1650#C  the  coefficient  of 
thermal  expansion  was  held  constant.  Figure  I I I - 1 
shows  the  variation  of  with  temperature  as  used 
in  this  study. 

c)  Yield  Stress  (cry) 

Available  data. for  yield  stress  for  HTS  ranged  from 

3 .4475  x  109  dynes/cm2  CS.S  x  104  psi)  to  4.3877  x 
9  2  4 

10  dynes/cm  (7.0  x  10  psi)  at  room  temperature. 

o  2 

The  low  value  of  3.4475  x  10  dynes/cm  was  used  for 
this  study.  Data  was  available  up  to  the  start  of 
the  transformation  range  at  S10°C.  The  curve  from 
this  point  on  was  estimated  based  upon  similar  data 
for  HY  steels  contained  in  References  6  and  7. 
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Since  no  dot*  for  hardening  modulus  (ET)  Was  available  -  ^ 

from  any  of  the  sources,  the  material  w»s  assumed  to 
be  elastic-perfectly  plastic.  However,  data  from  ] 

steels  with  similar  chemical  composition  and  heat  ] 

treatment  showed  significant  hardening  beyond  yield,  j 

at  least  at  room  temperature. 
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III-C.  Idealization  of  Material  Properties 


The  material  properties  used  in  the  NONSAP  analysis  are  tabulated 
in  Table  III-l.  NONSAP  is  limited  to  six  (6)  data  points  to 
define  material  properties  as  a  function  of  temperature.  The 
program  assumes  a  linear  relationship  between  the  data  points. 

Table  1 1 1 - 2  contains  the  material  properties  and  associated 
changes  in  temperature  used  in  the  Shrinkage  Force  Method  analy* 
sis.  The  material  properties  used  represent  the  average  for 
the  specified  ranges  in  Figure  I I I - 1 .  Based  on  the  shape  of 
the  curves  in  Figure  III-l,  three  average  values  for  a  and  E 
were  considered  a  reasonable  approximation.  Cool  down  analysis 
is  assumed  to  start  where  the  weld  metal  starts  to  achieve 
significant  strength  chosen  somewhat  arbitrarily  to  be  about 
1400*C. 
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IV.  APPLICATION  OF  THE  SHRINKAGE  FORCE  METHOD  (SFM)  TO 
WELDING  UNRESTRAINED  PLATES 


IV-A.  Introduction 

This  simplified  weld  distortion  prediction  technique  is  based 
upon  the  assumption  that  the  shrinkage  forces  developed  in  the 
weld  bead  during  thermal  cooldown  are  the  major  contributors 
to  distortion.  In  the  application  of  the  method,  the  geometry 
of  each  weld  bead  along  with  the  appropriate  temperature  changes 
(AT)  for  the  cooldown  of  the  bead  must  be  known.  In  its  simp¬ 
lest  form  the  linear  coefficient  of  expansion  (aT)  and  elastic 
modulus  (E)  are  used  to  calculate  the  linear  contraction  force 
bead  by  bead,  and  this  is  applied  to  the  workpiece  in  a  series 
of  numerical  solutions  for  distortions.  The  following  refine¬ 
ments  in  the  procedure  over  that  applied  in  Reference  2  were 
incorporated; 


a)  variable  (with  temperature)  material  properties  (E 
and  a  in  the  elastic  case) , 


b)  detailed  workpiece  and  weld  bead  geometry  repre¬ 
sentation, 

c)  actual  simulation  of  weld  bead  sequence, 


d) 
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remelt. 


The  advantages  of  the  mathod  include; 

a)  no  relatival/  expensive  time  dependent  nonlinear 
hea  ;  transfer  analysis  is  required, 

b)  not  restricted  to  the  simplified  conditions  of 
plain  strain  or  axisymmetric  behavior  in  the 
workpiece, 

c)  easy  and  economical  application  in  "real  world" 
shipyard  welding  situations. 

Although  a  heat  transfer  analysis  is  not  required,  it  may  be 
desirable  to  obtain  a  definition  of  the  final  solidus-liquidus 
boundary  as  a  result  of  each  weld  pass.  This  could  be  obtained 
from  a  heat  transfer  analysis,  empirical  methods  or  from 
experience.  Information  as  to  the  heat  affected  rone  is  also 
valuable,  but  not  required.  This  information  would  assist  the 
analyst  in  determining  the  volume  of  melted  material  available 
for  shrinkage  and  the  extent  to  which  the  high  temperatures 
have  altered  the  structure  a  d  thus  the  material  properties  of 
the  base  metal. 


i 
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IV-B.  Analytical  Approach 


The  various  stops  in  the  analysis  are  defined  in  the  flow  chart 
in  Figure  IV- 1.  The  flow  chart  is  general  in  that  it  includes 
the  steps  necessary  if  the  effects  of  remelt  are  considered  in 
the  analysis.  If  remelt  is  not  considered,  only  newly  deposited 
weld  metal  is  assumed  to  experience  thermal  shrinkage  during 
the  current  pass.  This  material  is  assumed  to  have  three  differ¬ 
ent  elastic  moduli  (E^,  i  «  1,2,3)  and  thermal  coefficients 
(aT^,  i  ■  1,2,3)  during  the  cooldown  that  occurs  in  the  three 
steps  (AT^,  i  -  1,2,3).  Values  for  B^,  aTi  and  AT^  are  given 
in  Section  III-C.  The  distortion  for  each  weld  bead  or  layer 
as  it  cools  down  is  computed.  When  there  is  no  account  for 
remelt,  these  displacements  are  simply  accumulated  and  added 
bead  by  bead  or  layer  by  layer.  The  model  reflects  the  changing 
stiffness  and  geometry  as  beads  or  layers  are  added.  Internal 
stresses  are  also  accumulated  and  added  throughout  the  process. 
When  remelt  is  considered,  the  internal  stresses  within  the 
previously  cooled  material,  which  is  then  assumed  to  reraelt,  are 
set  equal  to  zero  (i.e.  no  strength  is  assumed),  and  the  change 
in  deflection  is  computed  for  this  condition.  The  deflection 
and  internal  stresses  are  accumulated  for  solidification  and 
remelt  cycles  until  the  final  deformed  state  is  determined, 


Th«  approach  used  in  this  study  idealized  the  plat*  as  a  two 
dimensional  plana  strain  problem,  with  taro  strain  assuaad  in  tha 
diraction  of  walding.  As  discussad  in  Saction  II-C,  this  assump¬ 
tion  appaars  raasonabla  in  light  of  tha  tast  plata  gaoaatry  and 
axparimantal  daflaction  results.  All  models  employed  two 
diaansional  eight  noded  isoparametric  plane  strain  finite  ele¬ 
ments.  However,  the  method  is  not  restricted  to  this  geometry 
and  could  be  applied  to  axisymmetric  geometries  or  three  dimen¬ 
sional  solid  geometries  using  the  appropriate  finite  elements. 


IV-C.  Mathematical  Models 


Four  mathematical  models  were  developed  and  used  tp  predict  the 
distortion  caused  by  welding  the  first  side  of  the  unrestrained 
plate.  Model  1  consisted  of  a  six  layer  representation  of  the 
weld  metal.  Figure  IV- 2  shows  the  overall  model  of  the  plate 
and  the  detail  of  the  finite  element  model  in  the  vicinity  of 
the  weld  along  with  the  six  layers  of  weld  material.  One  half 
the  plate  and  weld  was  modeled  with  symmetrical  boundary  condi¬ 
tion  applied  along  the  line  B-F.  While  the  model  does  Include 
weld  joint  geometry,  it  does  not  consider  details  of  weld  bead 
geometry  or  remelt  effects.  From  previous  studies,  six  layers 
is  about  the  minimum  number  required  to  produce  "converged** 
final  displacements;  i.e.,  using  more  but  thinner  layers  would 
yield  essentially  the  same  results.  While  this  model  is  quite 
simplistic,  it  is  significant  because  it  reflects  the  types  of 
assumptions  that  for  most  practical  design  or  fabrication 
situations  one  might  be  forced  to  make;  i.e.,  the  details  of 
weld  bead  geometry  or  extent  of  remelt  are  not  normally  available. 
For  the  test  plate,  a  macroscopic  cross-section  of  the  weld  was 
photographed  at  S  1/2  power  magnification.  This  is  shown  in 
Figure  XV- 3.  With  this  data,  it  was  possible  to  roughly  esti¬ 
mate  the  weld  bead  geometry  and  volume  and  to  include  this  in 
formulating  model  2  which  is  shown  in  Figure  IV-4.  Due  to  remelt 
effects,  the  macroscopic  section  does  not  provide  a  direct  geo¬ 
metric  measure  of  the  'deposited  weld  metal.  Hence  a  good  deal 
of  judgement  an<l  estimating  were  required. 
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Model  2  does  not  assume  symmetric  response  about  the  centerline 
(£) .  Thus,  it  was  possible  to  study  tha  effect  of  non-symaetric 
application  of  wold  bands  on  tha  raaponsa.  Modal  3,  shown  in 
Figure  IV-SA,  considars  remelt  in  addition  to  wald  baad  geometry 
and  sequencing.  Tha  axtant  of  ramalt  was  estimated  from  tha 
macroscopic  saction.  Bands  1  and  2  ara  identified  in  Figura 
IV-  SB  with  tha  cross -hatchad  araa  idantifying  tha  ramalt  that 
occurs  in  Band  1  as  a  rasult  of  tha  daposlt  of  Baad  2.  Fig¬ 
ura  IV-6  shows  tha  antira  saquanca  of  weld  daposit  and  ramalt 
cyclas  usad  in  tha  analysis  to  simulate  tha  walding  of  this  joint. 

Modal  4  represents  a  change  in  tha  veld  groove  region  from 
Modal  3.  A  refined  wald  baad  geometry  was  usad  for  Model  4; 
this  is  shown  in  Figura  IV-7.  Tha  new  baad  geometry . is  based 
upon  rough  measurements  of  tha  fill  height  after  each  baad  was 
deposited.  Since  a  continuous  profile  of  upper  baad  surface 
was  not  taken  during  tha  test,  tha  baad  geometries  assumed 
represent  an  estimate  of  what  could  have  existed.  However,  this 
is  considered  an  improvement  over  tha  wald  baad  geometries  used 
in  Models  2  and  3.  A  second  change  between  Modal  4  and  Modal  3 
is  that  tha  cooldown  of  Baad  1  is  ignored,  and  the  analysis 
starts  with  tha  daposit  of  Baad  2.  Tha  experimental  results 
indicated  vary  little  distortion  for  tha  first  baad,  and  tha 
macrograph  (Figura  IV- 3)  provides  no  information  about  its  geo¬ 
metry,  Because  of  tha  lack  of  information  and  tha  small  amount 
of  distortion  associated  with  Bead  1,  neglecting  Bead  1  distortion 
appears  to  be  a  reasonable  approximation. 
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{  For  each  model,  the  finite  element  grid  in  the  weld  groove  for 

the  entire  first  side  weld  is  shown  in  the  figures.  In  fact, 
during  analysis,  the  model  includes  only  that  portion  of  the  weld 
which  exists  when  any  given  bead  or  layer  is  deposited  or 

s 

reaelted. 


IV" D.  Discussion  of  Results 


1)  Modal  1 

Figure  IV- 8  contains  plots  of  tip  def lection  versus  weld  bead 
number  or  corresponding  layer  number.  Also  shown  are  the  ex¬ 
perimental  results.  Overall  agreement  between  analytical  and 
experimental  results  is  reasonable  with  the  best  agreement 
occuring  for  the  first  beads  and  the  last  beads.  The  changes 
in  deflection  or  rotation  predicted  by  the  Intermediate  layers 
underestimates  the  distortion  produced  by  intermediate  beads  4 
and  S.  The  deflection  due  to  completed  first  side  weld  is 
underestimated  by  about  20%. 

< 

Figure  IV- 9  contains  a  plot  of  incremental  tip  deflections  due 
to  each  phase  of  the  cooldown  cycle  for  each  layer.  This  plot 
indicates  that  all  three  phases  of  the  cooldown  are  important, 
but  that  Phase  1  has  a  greater  influence  for  the  first  layers 
and  Phase  3  the  greater  influence  for  later  layers. 

Table  IV-l  compares  the  results  from  Reference  2  for  one  and 
two  constant  approximations  for  the  elastic  modulus  with  the 
results  from  Model  1  of  this  study  which  has  a  three  constant 
representation  of  modulus  with  temperature.  The  models  from 
Reference  2  utilised  a  one  constant  approximation  for  the 
coefficient  of  expansion  (a^.)  and  slightly  different  temperature 


changes  were  assumed.  The  table  does  show, -however,  the  improve* 
ment  of  the  distortion  results  using  the  simple  layer  approach 
with  more-  exact  representation  of  material  properties  with 
temperature.  The  final  experimental  deflection  for  welding  the 
first  side  is  also  contained  in  the  table. 

2)  Model  2 

Figure  IV- 10  contains  computed  incremental  and  cumulative  tip 
deflections  for  Model  2,  The  computed  incremental  deflections 
for  each  bead  show  best  correlation  with  experimental  data  for 
the  first  four  beads  and  the  last  two.  The  deflections  due  to 
beads  S,  6,  and  7  are  underestimated  by  approximately  SO  per 
cent.  The  cumulative  deflection  after  completion  of  the  first 
side  is  underestimated  by  approximately  30  per  cent,  i.e.  1.6 
inches  versus  1.1  inches. 

Computed  incremental  deflections  for  each  bead  and  each  phase  •  4 
of  cooldown  is  shown  in  Figure  IV- 11.  Consistent  with  the 
Model  1  results,  all  phases  of  cooldown  are  shown  to  be  impor¬ 
tant  contributors  to  deformation.  The  last  phase  of  cooldown 
(Phase  3)  has  the  most  significant  effect  for  all  beads  except 
number  three.  For  later  beads,  this  result  is  in  agreement 
with  results  for  Model  1.  However,  for  initial  beads  or  layers 
Model  1  showed  that  Phase  1  produced  the  largest,  deflections. 
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This  apparent  disagreement  batwaan  Modais  1  and  2  result s  is 
believed  to  ba  dua  to  differences  batwaan  tha  two  models  in 
geometric  representation  of,  in  particular,  tha  first  baad. 

Since  tha  plata  is  unrastraintd,  tha  daformation  will  ba  sansi- 
tiva  to  tha  geometry  assumed  for  its  first  baad  or  layar,  which 
is  quite  diffarant  for  Modais  1  and  2.  This  affact  should  ba 
considarably  lass  for  a  rastrainad  structura  and  in  fact  tha  nat 
affact  on  total  cumulative  daflactions  for  aven  this  unrestrained 
plata  is  not  particularly  significant.  Modal  2  rasults  showad 
no  significant  non- symmetric  response  of  tha  plata. 

3)  Modal  3 

Figure  IV- 12  contains  computed  incremental  and  cumulative  tip 
daflactions  for  Modal  3.  Tha  incremental  daflactions  for  each 
baad  agree  reasonably  well  with  experimental  data  except  for 
beads  1,  S  and  6.  Tha  computed  cumulative  deflection  for  com¬ 
pletion  of  tha  first  side  welding  agrees  within  a  few  par  cant 
of  tha  experimental  value,  i.a.  1,55  inches  versus  1.60  inches. 

Tha  Incremental  deflections  for  each  phase  of  the  cooldown  are 
given  in  Figure  IV-13.  Consistent  with  Model  1  and  2  .results, 
all  phases  are  relatively  important  with  Phase  3,  the  last 
phase  showing  somewhat  more  significance  for  Modal  3. 


The  incremental  deflections  due  to  remelt  ire  shown  in  Figure 
IV-14 .  These  results  show  that  remelt  causes  negative  deflec¬ 
tions  for  the  first  two  remelts  and  positive  deflections  for 
(  the  last  eight  remelts.  The  magnitude  of  these  effects  are 
small  in  relation  to  bead  cooldown  contributing  only  about  10 % 
to  final  deformation. 

4)  Model  4 

The  cumulative  tip  deflections  computed  with  Model  4  is  given  in 
Figure  IV-1S.  While  the  final  deflection  after  bead  9  does  not 
agree  with  experimental  value  quite  as  well  as  for  Model  3,  i.e. 

1.4  inches  versus  1.60  inches,  the  shape  of  the  curve  and  in¬ 
cremental  deflections  for  each  bead  are  judged  to  be  in  better 
overall  agreement.  This  overall  improvement  is  believed  to  be 
due  to  a  better  geometric  and  volume  approximation  to  each  bead. 
This  was  made  possible  by  use  of  both  the  macroscopic  section 
of  the  completed  weld  and  rough  measurements  of  the  depth  of 
fill  after  each  bead  was  laid  down.  As  stated  in  Section  IV-C, 
only  the  macroscopic  section  data  was  used  to  estimate  the  weld 
bead  geometry  and  volume  for  Models  2  and  3  which  is  not  sufficient 
for  an  accurate  representation. 

5)  Summary 

The  calculated  tip  deflection  for  Models  1  through  4  and  the 
measured  experimental  results  are  shown  on  Figure  IV-16.  Model  4 


is  considered  to  be  in  best  overall  agreement  with  experimental 
and  Model  2  the  worst.  However,  all  models  provide  reasonably 
good  estimates  of  the  final  deflection  for  the  completed  first 
side  weld. 


The  shape  of  the  deflection  curve  versus  bead  number  is  best 
approximated  by  Model  4  which  is  considered  to  have  the  most 
accurate  representation  of  actual  weld  bead  geometry,  and  which 
also  includes  the  effects  of  remelt. 

Model  1,  which  assumes  a  layered  weld  metal  deposit  and  no 
remelt  effects  ^stTl^tp'r^^des^reasonable  accuracy  for  this 

i  «  • 

problem. 

There  was  no  clear-cut  improvement  in  accuracy  obtained  in  going 
from  Model  1  to  Model  4,  even  though  each  model,  in  succession, 

.i 

appeared  to  be  better  than  the  previous  model.  Two  possible 
explanations  come  to  mind:  1)  the  detailed  buildup  of  the  weld 
bead  profile  is  not  known,  possibly  leading  to  an  erroneous 
conclusion  about  the  "goodness"  of  the  respective  models;  and 
2)  because  of  the  assumptions  and  simplifications  of  the 
"shrinkage  force"  analytical  approach,  as  applied  to  this  pro¬ 
blem,  the  type  of  model  refinement  attempted  here  may  not  be 
meaningful . 
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V.  APPLICATION  OF  THERMOPLASTIC  ANALYSIS  TO  WELDING 
UNRESTRAINED  PLATES 

V-A.  Introduction 

. — — — — in  ■■ 

Assuming  one  starts  with  a  validated  computer  code,  which 
possesses  the  necessary  analytical  capabilities  to  conduct  a 
welding* simulation  analysis,  the  accuracy  of  the  solution  will 
be  controlled  by  three  basic  inputs: 

1)  applied  loading 

2)  material  stress-strain  behavior  model 

3)  mathematical  discretization  of  structure. 

If  one  were  performing  a  static  linear  elastic  stress  analysis 
of  a  pressure  vessel,  the  first  two  inputs  are  straightforward 
and  the  only  real  work  involves  the  third  input.  However,  in 
welding-simulation  analysis  all  three  inputs  pose  real  problems 
for  the  analyst.  There  is  insufficient  experimental  data 
available  to  define  the  first  two  inputs  with  certainty.  The 
degree  of  experimental  sophistication  required  for  measurement 
of  welding- induced  phenomena  is  significantly  beyond  that 
required  to  measure  Young’s  Modulus.  Also,  because  of  the 
high  computer  run  cost  for  transient,  welding-simulation 
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analysts,  tht  number  of  dtgrtts  of  freedom  In  the  mathematical  f 

discretization  can  dictate  whether  an  analysis  is  attempted 
or  abandoned  as  too  costly.  The  effect  of  this  is  to  severely 
limit  the  conduct  of  modelling  studies  as  a  check  that  tho 
mathematical  discretization  of  the  structural  geometry  is 
adequate. 

« 

The  present  state-of-the-art  appears  to  be  1)  "best-guess"  time- 
dependent  temperature  distributions;  Z)  applied  to  a  structure 
with  "best-guess"  temperature-dependent  material  properties; 
and  3)  discretized  mathematically  by  an  economically  feasible 
number  of  degrees  of  freedom.  This  reality  makes  evaluation 
of  the  basic  analytical  capability,  as  represented  by  a  com¬ 
puterized  analytical  procedure,  extremely  difficult  since  lack 
of  agreement  with  actual  measured  response  could  be  caused  by 
incorrect , input ,  by  inadequate  theoretical  formulation,  or  by 
both.  Also  possible  is  the  situation  where  the  interaction 
of  a  number  of  inaccuracies  may  coincidentally  appear  to  produce 
accurate  results  when  compared  to  measured  data. 

With  these  considerations  in  mind,  application  of  the  NONSAP 
computer  code  to  transient,  thermo -plastic  analysis  of  welding- 
induced  behavior  should  be  viewed  more  as  a  feasibility  and 
qualitative  evaluation,  than  as  a  quantitative  evaluation  for 
correlation  to  experimental  results. 

i 
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The  NONSAP  finite  element  program  was  modified  in  1977  to  per- 
form  thermoplastic  analysis  with  temperature  dependent  material 
properties  (Reference  2) .  Additional  modifications  were  made 
in  1978  to  assist  in  the  analysis  of  weld  distortion  problems. 
Specifically,  a  "birth"  and  "death*'  option  has  been  added  for 
the  two  dimensional  finite  elements  representing  the  weld 
deposit.  At  predetermined  times  in  the  analysis,  elements  in 
the  weld  groove  region  can  be  born  (achieve  strength  from 
cooldown)  or  die  (melt).  This  feature  assists  in  modelling 
of  the  weld  bead. 

One  of  the  difficulties  encountered  in  the  analysis  of  unrestrained 
butt  welded  plates,  supported  on  knife  edges,  is  that  the  plates 
are  basically  free  to  rotate  as  rigid  bodies  during  the  weld 
process.  Thus,  during  welding,  the  plates  lift  off  one  of 
the  2  sets  of  supports;  therefore,  a  nonlinear  support  with 
large  compressive  stiffness  and  essentially  zero  tensile  stiff¬ 
ness  had  to  be  simulated.  This  was  accomplished  using  1-D 
truss  finite  elements  with  elastic,  nonlinear  material  properties. 
With  the  required  program  development  completed,  the  analysis 
of  the  experiment  described  in  Section  III  was  initiated. 

Transient  temperature  distributions  in  the  plates  for  deposit 
of  weld  beads  1  and  2  were  calculated  in  1977  using  the  TEMP 
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computer  coda  (saa  Reference  2).  These  vara  used  in  1978  as 
input  to  NONSAP  for  tha  analysis  of  distortion  causad  by  daposit 
of  assuaad  wald  layers  1  and  2. 

Tha  wald  layar  geometry  for  the  NONSAP  analysis  is  tha  same 

as  that  used  in  GENSAM  Modal  1,  described  in  tha  previous 

) 

section.  A  detailed  weld  bead  definition,  such  as  was  used 
in  GENSAM  Models  2,  3,  and  4,  was  not  attempted  in  tha  NONSAP 
analysis.  The  lack  of  detailed  weld  bead  description  in  the 
NONSAP  model  increases  the  difficulty  of  interpreting  and  eval¬ 
uating  the  results  obtained,  since  there  is  not  a  one-to-one 
correspondence  between  the  analytical  model  and  the  experimental 
model  for  either  weld  layer  1  or  weld  layer  2.  This  choice 
of  weld  bead  description  was  made  based  on  concern  for  com¬ 
puter  run  cost  and  the  limitation  on  degrees  of  freedom  in 
NONSAP,  which  uses  an  in-core  solution  technique.  More  detailed 
(and  correspondingly  more  expensive)  models  could  be  analysed 
using  the  ADINA  computer  code,  which  is  a  general  purpose 
code,  based  on  NONSAP,  developed  by  Prof.  Bathe  at  MIT. 


V-B  Analysis  Methods 


Figure  Y-l  shows  the  overall  NONSAP  finite  element  model  used 
for  the  thermoplastic  analysis.  The  simple  supports  at  A 
and  B  have  been  replaced  with  a  spring  trusswork.  These  sup- 
ports  have  the  ability  to  roll  freely  in  the  Y  direction  and 
constrain  translation  in  the  -2  direction.  This  is  accomplished 
by  specifying  nonlinear  properties  for  the  truss  members. 

This  mechanism  for  the  supports  generally  worked  well.  Fig- 
ure  V-2  shows  the  model  of  the  weld  groove  region  used  In  the 
NONSAP  analysis.  Translational  motion  in  the  Y  direction  was 
constrained  at  the  weld  centerline  (AB) .  Both  Layers  1  and  2 
are  shown  in  Figure  V-2;  however,  their  births  occur  at 
different  times  during  the  analysis  and  their  mutual  presence 
in  Figure  V-2  is  purely  illustrative. 

Four-noded  isoparametric,  quadrilateral,  plane-strain  finite 
elements,  with  a  2x2  Gaussian  integration  scheme,  were  util¬ 
ised  in  the  analysis.  The  elements  were  reduced  to  triangular 
shapes  where  necessary,  to  represent  the  geometry  and  for 
transition  from  large  element  sUo  to  small  element  site  as 
the  weld  region  is  approached. 
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While  the  grid  is  significantly  refined  in  the  weld  region, 
it  is  still  too  crude  to  represent  the  detailed  thermo-plastic 
stress  behavior  of  the  weld  metal  and  adjacent  plate.  To 
accurately  predict  these  stresses  during  and  after  welding, 

the  element  density  would  have  to  be  5-10  times  greater  for 

\ 

the  elements  used.  However,  as  a  demonstration  problem,  to 
determine  the  overall  distortion  behavior  of  the  plate,  the* 
element  density  was  judged  adequate. 

The  only  load  on  the  plate,  other  than  the  transient  temperature 
distribution  caused  by  welding,  is  the  dead  weight  of  the 
plate  itself. 

From  the  transient  heat  transfer  analysis,  a  file  of  nodal 
temperature  distributions  as  a  function  of  time  were  generated 
to  simulate  the  heatup-cooldown  cycle  caused  by  deposition  of 
weld  layers  1  and  Z.  For  the  NONSAP  analysis,  a  sequence  of  : 

‘l 

temperature  distributions  was  chosen  from  this  output,  such 
that  the  maximum  temperature  change  at  any  node  in  going  from 
time  (l)  to  time  (1*1)  was  j AT j  <_  2G0*C. 

From  the  nodal  temperatures,  the  temperature  at  the  four  Gauss* 
points  in  each  finite  element  is  calculated  in  NONSAP.  The 


tangent  stiffness  is. then  calculated,  using  the  material  pro¬ 
perties  associated  with  the  current  temperatures  and  stresses 
at  the  Gauss  points.  The  increment  of  Thermal  Load  in  going 
from  time  (i)  to  time  (i+1)  is  applied  to  obtain  the  increment 
in  nodal  deflections  and  the  associated  incremental  stresses 
at  the  Gauss  points,  assuming  elastic  behavior.  Then  the  total 
stress  at  each  Gauss  point  is  compared  to  the  Von  Mises  yield 
surface  at  the  current  temperature  to  determine  whether  the 
increment  in  stress  calculated  elastically  is  admissible  or 
whether  the  yield  surface  has  been  exceeded.  If  the  yield 
surface  has  been  exceeded,  the  stresses  are  adjusted  to  be 
consistent  with  the  specified  plastic  behavior  model. 

Equilibrium  iterations  are  performed  until  the  displacement 
solution  has  converged  to  within  a  .005  error  tolerance.  Within 
each  iteration,  the  stresses  are  updated,  checked  against  the 
yield  surface,  and  corrected  if  necessary. 
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V-C.  Analysis  Results 


Figure  V-S(a)  shows  the  average  temperatures  o£  Weld  Layers  1 
and  2  during  the  analysis  and  Figure  V-3(b)  shows  the  tip 
deflection  of  the  plate  during  the  analysis.  The  solution  for 
Layer  1  deposit  had  25  load  (temperature)  steps  and  the  solution 
for  Layer  2  deposit  had  44  load  steps. 

The  initial  approach  was  to  analyze  one  temperature  sequence, 
including  the  two  heatup-cooldown  cycles  and  introducing  the 
Weld  Layers  1  and  2  at  the  appropriate  times.  This  approach 
had  to  be  abandoned  when,  during  the  heatup  for  Layer  2,  un- 
resolvable  convergence  difficulties  developed.  As  the  temper¬ 
ature  of  the  previously  deposited  Layer  1  increased,  during 
the  2nd  heat-up,  the  yield  strength  of  Weld  Layer  1  decreased 
to  the  point  where  the  plate  deflection  could  no  longer  be 
supported,  and  consequently  the  deflection  rapidly  reduced 
toward  zero.  A  displacement  solution  showing  the  plate  settled 
back  on  the  supports  was  obtained,  but  the  reaction  forces  at 
the  supports  were  grossly  in  error.  While  considerable  effort 
was  made  to  obtain  convergence  of  the  support  reactions,  no 
solution  startegy  was  found  to  accomplish  this.  When  the 
solution  was  continued  to  the  next  load  3tep,  the  displacement 
solution  blew  up. 
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At  this  point,  it  was  decided  to  start  a  new  analysis  at  the 
beginning  of  the  second  heatup-cooldown  cycle  for  Layer  2. 

This  approach  is  not  rigorous  in  that  the  stresses  are  assumed 
to  be  zero  at  the  beginning  of  the  2nd  analysis,  while  in 
reality  there  is  a  non-zero  stress  field  resulting  from  the 
Layer  1  heatup -cooldown  cycle.  However,  because  these  stresses 
in  Layer  1  are  relieved  during  the  second  heat-up  cycle  as 
Layer  1  is  effectively  re-melted,  this  shortcoming  does  not 
appear  to  be  particularly  significant.  Indeed,  as  can  be  seen 
from  Figure  V-3(b),  Analysis  2  showed  that  an  initial  tip 
deflection  developed,  but  could  not  be  maintained  as  the  tem¬ 
perature  increased;  and  the  plate  settled  back  on  the  supports 
at  about  the  same  point  as  in  Analysis  1.  Error  in  the  support 
reactions  also  occurred  at  this  point  of  Analysis  2;  on  the 
next  load  step,  howover,  balance  between  the  weight  of  the 
plate  and  the  support  reactions  was  re-established. 
ii 

From  Figure  V-3(a),  it  can  be  seen  that  Weld  Layer  1  is  effec¬ 
tively  re-melted  during  the  heatup  for  the  deposit  of  Weld 
Layer  2  -  achieving  an  average  temperature  of  17S0°C  at  the 
point  where  Weld  Layer  2  is  introduced.  As  noted  above,  this 
results  in  the  plate  deflection,  returning  to  zero  before  Weld 
Layer  2  is  deposited.  Effectively,  then,  the  analysis  of  Weld 
Layer  2  approaches  the  case  of  the  deposit  of  an  initial  weld 
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layer  whose  size  is  the  sum  of  Weld  Layers  1  and  2.  This 
phenomenon  may  explain  whey  the  final  predicted  tip  deflection 
Is  almost  the  same  after  the  cooldown  of  Layer  2  as  if  was 
after  the  cooldown  of  Layer  1. 

On  cooldown,  once  sufficient  strength  has  been  achieved  by 
the  wild  metal  to  provide  the  net  moment  necessary  to  deflect 
the  plnte,  this  net  moment  remains  constant  and  the  lateral 
deflection  of  the  plate  is  controlled  by  the  thermally- induced 
contraction.  As  a  check  on  the  NONSAP  solutions,  the  moment 
that  must  exist  at  the  inside  support,  to  allow  lift-off  from 
the  outside  support,  was  calculated  and  compared  to  the  net 
moment  duo  to  the  axial  Jtre«s  distribution  predicted  by  NONSAP 
at  the  inside  support.  The  moment  due  to  predicted  stresses 
was  approximately  equal  to  the  calculated  required  moment  for 
all  load  steps  where  the  plate  is  deflected,  and  was  less 
than  the  required  moment  for  all  load  steps  where  the  plate 
is  undeflected. 

The  plate  behavior  for  deposit  of  Weld  Layer  1  is  straight¬ 
forward.  Until  the  finite  element  representing  Layer  1  has 
cooled  to  a  temperature  low  enough  to  achieve  sufficient  strength, 
the  plate  remains  undeflected.  At  760*C,  the  strength  of  the 
metal  starts  to  increase  dramatically  (see  Figuro  III-l),  and 


as  th«  mstal  cools  further,  lift  off  frost  the  outside  support 
begins  (Load  Step  21).  From  this  point  (T*S70*C) ,  the  tem¬ 
perature  in  and  adjacent  to  Weld  Layer  l  decreases  uniformly 
until  room  temperature  is  reached.  The  thermal  contraction 
of  Weld  Layer  1  and  adjacent  base  plate  controls  the  plate 
deflection.  Because  the  weld  layer  lies  above  the  neutral 
axis  of  the  plate,  the  plate  rotates  such  that  the  tip  deflects 
upward.  The  final  predicted  tip  deflection  was  0.116  cm. 

The  plate  behavior  for  deposit  of  Weld  Layer  2  is  more  complex. 
From  Figure  V»3(b),  Analysis  2  shows  three  regions  during  the 
heatup- cooldown  cycle  where  the  plate  deflected.  In  Analysis  2, 
Weld  Layer  1  is  in  place  and  is  initially  subjected  to  an. 
increase  in  temperature.  A  uniform  temperature  increase  of 
Weld  Layer  1  would  be  expected  to  produce  a  negative  tip 
deflection.  However,  the  NONSAP  analysis  predicted  a  positive 
tip  deflection  during  initial  heatup.  To  explain  this 
behavior,  it  is  necessary  to  examine  the  temperature  distri¬ 
bution  in  the  element  representing  Weld  Layer  1  and  the  adjacent 
element  in  the  base  plate. 

The  transient  heat  transfer  analysis  predicted  temperatures 
at  the  locations  of  the  nodes  in  the  NONSAP  model.  The  temper¬ 
atures  at  the  Gauss  pts.  in  each  element  are  a  linear  combination 
of  these  nodal  temperatures.  Figure  V-4  shows  the  predicted 
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nodal  temperatures  and  calculated  temperatures  at  the  Gauss 

> 

pts.  for  the  element  representing  Weld  Layer  1  and  the  adjacent 
element  in  the  base  plate  at  load  step  13  (where  the  early 
plate  deflection  peaks). 

The  nodal  temperatures  show  a  sharp  thermal  gradient  from  the 
weld  centerline  into  the  base  plate.  At  the  top  of  Weld 
Layer  1  and  the  adjacent  element,  the  three  nodal  temperatures 
approximate  the  highly  nonlinear  thermal  gradient.  At  the 
bottom  of  Weld  Layer  1  and  adjacent  element,  the  two  nodal 

temperatures  can  only  represent  a  linear  thermal  gradient 

* 

over  the  same  distance  into  the  plate. 

The  temperatures  at  the  Gauss  pts.  in  these  two  elements  are 
calculated  by  linear  interpolation  from  the  nodal  temperatures. 
Examination  of  the  Gauss  pt.  temperatures  indicates  that  the 

restrained  thermal  growth  will  be  greater  at  the  bottom  of 

(1 

Wold  Layer  1  and  adjacent  element  than  at  the  top.  As  a  re¬ 
sult,  the  plate  is  rotated  such  that  the  tip  deflects  upward. 

Based  on  this  evaluation  of  the  plate  behavior  during  early 
heatup,  it  is  concluded  that  the  positive  plate  deflection 
is  not  caused  by  physical  considerations,  but  rather  by  the 
mathematical  modelling  approximations  to  the  physical  problem. 
The  finite-elemont  grid  is  inadequate  to  represent  the  actual 
thermal  gradient  over  this  short  distance. 
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This  doflection  is  transient;  it  is  totally  reversed  when  Weld 
Layer  1  loses  the  strength  to  support  the  plate  in  a  deflected 
position,  as  the  temperature  increases  above  «6G0*C.  There¬ 
fore,  this  early  plate  deflection  does  not  contribute  to  the 
final  permanent  plate  deflection  after  cooldown  of  Weld  Layer  2. 

The  discussion  above  also  applies  to  the  plate  deflection 
behavior  predicted  in  Analysis  1,  during  the  second  heatup 
stage  (dashed  line  in  Figure,  V- 3(b)) .  Here,  the  additional 
positive  tip  deflection  is  much  less,  because  stresses  approach¬ 
ing  yield  stress  already  exist  after  the  cooldown  of  Layer  1. 

The  thermally  induced  doflection  behavior  is  not  stress- free, 

« 

sines  the  temperature  increase  is  far  from  uniform.  The  cap¬ 
ability  to  carry  additional  stress  is  limited  by  the  pre-existing 
stresses,  which  limits  the  amount  of  additional  positive  tip 
deflection. 

ii 

Further  along  in  Analysis  2,  after  settling  back  on  the  sup¬ 
ports,  the  plate  remains  undeflected  through  the  remaining 
heatup,  the  birth  of  Wvld  Layer  2,  and  the  initial  cooldown. 

Above  1SOO*C,  the  metal  has  negligible  strength  and  cannot 
support  a  plate  deflection.  At  load  step  31,  when  the  average 
temperature  of  Weld  Layers  1  and  2  drops  below  iS00*C,  the 
plate  tip  starts  to  deflect  [see  Figure  V-3(b)].  Because  the 
total  thickness  of  the  weld  buildup  is  now  doubled,  £he 
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moment  needed  to  deflect  the  plate  can  develop  at  a  much  lower 
material  strength  level ,  and  a  correspondingly  higher  temper¬ 
ature.  The  transient  plate  deflection  between  load  steps  31 
and  39  appears  to  be  caused  by  the  different  cooling  rates 
for  Wold  Layers  1  and  2  [see  Figure  V-3(a)].  Initially,  Weld 
Layer  2  cools  at  a  faster  rate  and  has  a  lower  avorage  tem¬ 
perature  than  Weld  Layer  1,  causing' an  upward  tip  deflection. 
Then,  Weld  Layer  1  cools  at  a  faster  rate  than  Weld  Layer  2 
and  the  average  temperatures  converge, causing  a  downward  tip 
deflection.  Finally,  by  load  step  39,  the  cooling  rates  and 
average  layer  temperature  equalise. 

This  transient  deflection  occurs  in  the  weld  layer  temperature 
range  1500"0760<,C,  before  the  major  increase  in  material 
strength  occurs  (see  Figure  X 1 1 - 1 ) .  While  the  deflection  is 
plausible,  based  on  the  weld  layers  temperature  behavior, 
its  physical  significance  is  difficult  to  assess.  Transient 
experimental  data  is  unavailable  for  correlation;  and  it  is 
conceivable  that  lack  of  restraint  on  the  plate,  coupled  with 
approximations  in  the  heat  transfer  and  structural  models, 
might  produce  an  artificial  effect  in  the  high-temperature 
region  of  the  solution.  However,  whether  physical  or  not, 
this  transient  deflection  does  not  contribute  to  the  final 
plate  deflection  after  complete  cooldown. 
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The  development  of  the  final  plate  tip  deflection  begins  with 
load  step  40,  when  the  weld  layer  temperatures  drop  below  760°C. 
From  load  steps  40  (T*670°C)  to  44  (room  temperature),  the 
temperature  decreases  uniformly  in  and  adjacent  to  Weld  Layers  1 
and  2.  The  pattern  of  deflection  growth  is  very  similar  to 
that  predicted  in  Analysis  1  for  cooldown  of  Wold  Layer  1. 

The  final  tip  deflection  after  cooldown  of  Weld  Layer  2  is 
0.117  cm.  upward,  which  is  about  II  higher  than  the  deflection 
after  cooldown  of  Weld  Layer  1. 

It  is  difficult  to  evaluate  the  accuracy  of  the  final  tip 
deflection  predictions  from  Analyses  1  and  2.  Direct  compar¬ 
ison  to  experimental  data  is  not  possible  because  there  is  no 
known  correspondence  between  assumed  Weld  Layers  1  and  2  and 
the  actual  weld  beads  from  the  experiment.  However,  limited 
data  was  taken  during  the  experiment,  showing  the  average 
level  of  joint  fill  after  beads  3,  4,  6,  8,  and  9.  Unfortunately , 
similar  measurements  were  not  taken  after  beads  1  and  2. 

Figure  V'-S  shows  a  plot  of  experimental  tip  dofloction  vs. 
height  of  joint  fill,  extrapolated  back  to  zero  deflection 
at  zero  fill  height.  Also  shown  are  the  NONSAP  results  for 
assumed  Weld  Layers  1  and  2.  While  the  NONSAP  results  do  not 
show  an  increase  in  tip  deflection  from  Weld  Layer  1  to  Weld 
Layer  2,  the  predicted  results  are  reasonably  close  to  the 
extrapolated  experimental  curve. 
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Without  additional  heat  transfer  and  tharmo- structural  analysis 
results  for  tha  remaining  veld  layers,  the  NONSAP  results 
reported  herein  must  be  considered  inconclusive.  Future  appli¬ 
cations  of  the  NONSAP  approach  should  consider  the  complete 
fill  of  the  weld  joint,  to  better  assess  the  accuracy  of  the 
analytical  predictions. 
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FIGURE  31-5  NONSAP  v*.  EXPERIMENTAL  DEFLECTIONS 


VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


VI-A.  Shrinkage  Force  Analysis 

The  developments  incorporated  in  this  approach  this  year 
have  brought  about  an  overall  improvement  in  the  ability  to 
predict  distortion.  A  simple  layer  by  layer  approach  with¬ 
out  remelt  effects  provides  a  good  initial  solution  for  the 
distortion  problem.  Although  the  effects  of  remelt  were  not 
included  in  model  1  (layer  model),  results  with  later  models 
indicated  that  remelt  reduced  the  distortion  in  the  early 
beads  and  becamo  additive  in  the  later  beads.  Applying  this 
effect  to  the  results  for  model  1  in  Figure  IV- 8,  it  would 
improve  the  overall  prediction  of  model  1.  The  bead  approach 
used  in  models  2,  3,  and  4  does  provide  a  better  approximation 
to  the  distortion.  However,  the  information  required  to  im¬ 
plement  these  models  is  extensive  (bead  profiles,  remelt 
regions)  and  there  remains  a  problem  employing  models  of  this 
type  on  shipyard  probloms. 

For  the  unrestrained  plate  problem,  reasonably  accurate  estimates 
of  distortion  were  obtained  without  considering  plasticity 
effects  in  the  weld  metal.  This  is  not  surprising  since  the 
structures  being  welded  together  offer  minimal  restraint  to 
the  shrinking  weld  metal.  However,  it  is  not  clear  why  the 
predicted  distortions  were  less  than  the  measured  distortions, 
since  the  assumption  of  elastic  behavior  of  the  weld  metal 
would  be  expected  to  yield  an  upper  bound  analytical  prediction 
for  distortion. 
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The  greatest  discrepancy  in  the  results  was  the  inability 
of  any  of  the  models  to  predict  the  general  increase  in  deflec¬ 
tion  that  the  experiment  showed  in  going  from  Bead  4  to  Bead  6 
(Figure  IV- 16).  No  explanation  can  be  found  for  this  except 

that  the  bead  representation  and  individual  bead  remelt  geo¬ 
metry  are  still  not  sufficiently  defined.  Because  of  the 
required  detailed  bead  description,  application  of  models  2, 

3  and  4  is  somewhat  limited.  Methods  which  provide  for  pre¬ 
dicting  bead  geometry  and  remelt  will  have  to  be  developed 
before  models  such  as  these  can  be  utilized  fully.  This 

appears  to  be  feasible  if  the  parameters  of  the  welding  pro- 

* 

cess,  such  as  rate  and  volume  of  deposit,  average  penetration 
etc.,  are  known. 

Based  on  the  results  of  the  SFM  study  the  following  conclusion 
can  be  drawn: 

a)  the  simple  layer  model  for  unrestrained  welds 
provide  an  adequate  prediction  of  distortion  for 
shipyard  problems. 

b)  the  effects  of  material  property  variation  with 
temperature  are  important  even  in  the  simplest 
model  of  weld  distortion. 


c)  plasticity  effects  were  not  significant  in  predicting 
weld  distortion  of  unrestrained  plates;  however 
plasticity  affects  In  restrained  welds  will  be  very 
important  for  accurate  prediction  of  weld  distortion. 

d)  the  SFM  retains  several  advantages  (see  Section  IV-A) 
over  more  sophisticated  methods,  and  is  not  re* 
stricted  to  a  particular  welding  process  or  material. 

e)  useful  results  were  obtained  by  SFM  for  this  problem 
even  with  a  coarse  finite  element  subdivision  of  the  weld 

region;  however,  the  analysis  of  restrained  welds 

* 

with  plasticity  effects  may  require  finer  grids  to 
insure  sufficiently  accurate  calculation  of  the 
tangent  stiffness  in  the  weld  region. 
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VI "B.  NONSAP  Analysis 


Based  on  the  NONSAP  results  obtained  in  thif  study,  the  per* 
manent  distortion  appears  to  be  controlled  by  behavior  below 
760*C,  where  the  temperature  decrease  in  and  adjacent  to  the 
deposited  weld  metal  is  essentially  uniform.  Transient  heat 
transfer  analysis  is  important  to  define  the  extent  of  melting 
of  solid  metal  when  the  current  weld  bead  is  being  deposited. 
However,  if  this  could  be  reasonably  estimated  based  on  past 
experience  and  experimental  studies,  the  heat  transfer  analysis 
may  not  be  essential  for  distortion  prediction.  Also,  mater¬ 
ial  property  variations  with  temperature  above  760*C  -  the 
region  where  they  are  most  difficult  to  define  -  may  not  be 
essential  to  the  analysis  of  distortion.  These  concepts  are 
currently  being  investigated  under  an  Electric  Boat  Division 
IRAD  project  to  develop  a  second  generation,  shrinkage  force 
method  for  predicting  weld  distortion. 

As  discussed  in  Section  V,  the  NONSAP  predictions  of  distortion 
obtained  for  assumed  Weld  Layers  1  and  2  are  difficult  to 
evaluate  against  the  available  experimental  data;  the  results 
are  considered  inconclusive  from  a  quantitative  viewpoint. 
However,  conceptually,  the  results  are  promising  because  the 
analytical  model  has  demonstrated  the  ability  to  simulate 
fundamental  behavior  occurring  during  the  welding  process 
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With  improved  input,  refined  modelling,  and  more  comprehensive 
comparison  with  experimental  data,  the  NONSAP  approach  may 
prove  to  be  a  valuable  analytical  tool  for  prediction  of  welding- 
induced  structural  behavior. 

Considerable  work  remains  to  be  done  in  both  the  heat  transfer 
and  structural  analysis  areas  before  detailed  residual  stresses 
can  be  accurately  predicted.  Based  on  the  knowledge  gained 
in  the  present  study,  the  following  conclusions  are  drawn  con¬ 
cerning  residual  stress  prediction: 

a)  performance  of  transient  heat  transfer  and  thermo¬ 
plastic  analyses  for  multi-pass  welding  problems 
is  expensive,  and  for  most  shipyard  applications 
it  would  be  difficult  to  economically  justify  this 
approach  for  distortion  prediction;  for  detailed 
residual  stress  analysis,  this  approach  not  only 
appears  to  bo  necessary  but  also  will  probably 
cost  an  order  of  magnitude  more  than  a  distortion 
analysis  for  comparable  accuracy. 

b)  nonlinear  heat  transfer  analysis  and  the  correct 
determination  of  the  transient  temperature  solutions 
is  essential  to  welding  residual  stress  analysis; 
more  attention  should  be  directed  toward  improvements 
in  this  area. 
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accurate  representation  of  material  property  varia¬ 
tions  with  temperatures  up  to  melt  appears  necessary 
for  residual  stress  analysis;  additional  work  is 
needed  in  this  area. 

to  accurately  predict  temperatures  and  stresses  in 
the  vicinity  of  the  weld  groove,  the  finite  element 
discretization  must  be  fine  enough  to  pick  up  sig¬ 
nificant  gradients  over  distances  on  the  order  of 


cm. 


VI -C.  Recommendations 


The  following  tasks  are  recommended  for  future  investigations: 

1)  Analytically  investigate  the  effects  of  varying 
key  parameters  in  the  heat  transfer  problem  (heat 
deposit  mechanism,  heat  dissipation  to  environment, 
material  properties  as  functions  of  temperature) 

on  the  predicted  transient  temperature  distributions 

» 

that  occur  during  girth  welding  of  HY-80/130  steel 
cylinders. 

2)  Analytically  predict*  transient  temperature  solution 
for  girth  welding  of  HY-130/HY-80  steel  cylinders. 

3)  Compare  analytical  temperature  predictions  from  2) 

with  the  experimental  results  obtained  by  Professor 
Masubuchi  at  MIT.  " 

4)  Analytically  predict  transient  thermoplastic  distort 
tions  for  the  girth  welding  of  HY-130/HY-80  cylinders 
and  compare  with  experimental  results  from  MIT. 

5)  Apply  the  Shrinkage  Force  Method  (SFM)  to  predictions 
of  thermomechanical  distortions  resulting  from  girth 
welding  of  HY-150/HY-80  cylinders  and  compare  with 
experimental  results  from  MIT. 
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6)  Apply  the  Shrinkage  Force  Method  to  prediction  of 
thermomechanical  distortions  resulting  from  girth 
welding  of  HY-80  submarine  cylinders  under  shipyard 
conditions.  Compare  this  data  with  measurements 
taken  during  actual  weld  construction. 

7)  Develop  a  detailed  pXan  for  determination  of  mattsrial 
property  data  for  high  strength  steels  such  as 
HY-80. 


as 
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